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Spin transport properties of the one-dimensional Heisenberg antiferromagnetic spin systems for 
botli 5* = 1/2 and 5* = 1 are studied by applying twisted boundary magnetic field. The spin current 
displays significantly different behavior of the spin transport properties between 5* = 1/2 and 5 = 1 
cases. For the spin-half case, a London equation for the current and the detection of an alternating 
electric field are proposed for the linear response regime. The correlation functions reveal the spiral 
nature of spin configuration for both ground state and the spinon excitations. For the spin-one 
chain otherwise, a kink is generated in the ground state for the size is larger than the correlation 
length, leading to an exponential dependence of spin current with respect to the chains length. The 
midgap state emerges from the degenerate ground state even for small boundary fields. 
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A significant amount of experimental and theoretical 
efforts has been focused on the controlling of spin degree 
of freedom in recent years. P, 0,0 For the most part the 
research effort has been concentrated on the dilute mag- 
netic semiconductors as well as the spin Hall effect 
emergingfrom the spin-orbital coupling in two dimen- 
sions 0,MIH0- However, spin transport properties in 
pure spin systems, such as the Heisenberg system, is also 
of great interests because many of the novel concepts 
associated with spin conduction can be tested without 
the interference of charge degrees of freedom |^ 0, llo| . 
For instance, ballistic transport characterized by a finite 
Drude weight or spin stiffness has been found for both 
the integrable systems and certain class of Luttinger 
liquids. y Recently, Meier and Loss studied the magne- 
tization transport properties in a finite spin-half Heisen- 
berg chain linked to two bulk magnets. [ij They obtained 
a finite spin conductivity for a confined antiferromag- 
netic chain and predicted that a magnetization current 
produces an electric field. Alternatively, Schiitz et al[l^ 
investigated a mesoscopic spin ring in the inhomogcncous 
magnetic field to search for persistent spin current for dif- 
ferent spins. On the experimental front, mean free paths 
of several hundred nanometers were found to suggest a 
quasi-ballistic transport of one dimensional elementary 
spin excitations in Sr2Cu03 and SrCu02 samples 0. 

In the present work, we study the spin transport 
properties of one-dimensional (ID) Heisenberg antifer- 
romagnetic (HAF) spin models. It is well known that 
the integer-spin HAF chain is distinguished by a finite 
gap in spectrum from the half- integer spin one. [3 As 
expected that the spin transport is diffusive for the 
integer spin case, experimental NMR^^ and thermal 
conductivity |0| measurements have indicated finite spin 



diffusion and thermal diffusion constants in AgVP2S6. 
On the other hand, Fujimoto based on the integrability 
of the nonlinear cr-model suggested that the spin trans- 
port is ballistic in the perfect ID spin one system. 
Therefore, a general criteria on ballistic spin transport 
are yet to be established generally for either homogenous 
spin chains with different spins or mesoscopic (quasi-)one 
dimensional Heisenberg systems with an inhomogeneous 
magnetic field. This motivates us to devise a simplified 
case in which twisted magnetic fields are applied to both 
ends of a finite HAF chain with both 5 = 1/2 and S = \. 
Whether the twisted boundary magnetic fields, breaking 
the translational as well as SU(2) symmetry, can prop- 
erly drives the spin to flow through the chains alterna- 
tively presents distinguishable transport nature between 
the spin half and one cases. 

The Hamiltonian we will consider reads: 



N-l 



H — J Si ■ S,;+i — hi • Si — h2 • Spf, 



(1) 



where is the spin operator for cither S = 1/2 or 1 at 
the ith site, respectively, hi and hjv are the magnetic 
fields applied to the spins Si and Sat. For convenience, 
we set J = 1 and |hi| = |hjv| = h unless specified, take 
hi = (0, 0, h) and hAr = {h sm9, 0, h cos6) where 9 is the 
angle between hi and hjv and in [0, tt]. We note that such 
a system can be realized experimentally by attaching two 
magnetic leads on two ends of the chains. [ll] 

We employ the density matrix rcnormalization group 
(DMRG) methodilX, 18,19] to study the present twisted 
magnetic field effects. In our computations, wc use both 
the infinite and the finite size algorithm. The maximum 
number of sweeps is five and the number of states is kept 
up to 400. The efficiency of computations is not reduced 
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significantly by the absence of tlie conserved S**"*"'. Tlie 
truncation errors are about 10~^^ and the relative errors 
maintained below one percent as examined by increas- 
ing more kept states and sweeps. The computations are 
performed up to 100 sites usually and 400 sites for the 
proper extrapolations needed for thermodynamics limit. 

It is natural to first examine the twisted magnetic field 
effects on the correlation function 



(2) 



of 6* at /i = 1. For 5=1/2, we found that 



CUr) = {-ir{SlSU. 



for the spin z— component and its alternative form 

GUr) = C.M~{-ir{Sl){SU,). (3) 

Figure shows the correlation function Czz {r) with re- 
spect to the distance r for the chain length L = 100 and 
various angles at /i = 1. We found that Czz{r) is nonzero 




FIG. 1: Spin correlation function Czzij') for both 5=1/2 
(left penal) and 5 = 1 (right penal) for six angles 9 = Q, 7r/6, 
7r/3, 7r/2, 27r/3 and tt with h = l. The insets show Gzz{r). 

almost for all r and 6 for the S = 1/2, while Czzir) is 
finite only when r is nearby both ends of the chain owing 
to gapful excitation for the S = 1 case. This is intrin- 
sically associated with their algebraic and the exponent 
behaviors at /i = for S" = 1/2 and 1, respectively. As 
expected, C^^(r, 9) + Czz{r, -k ~ 9) = 2Czz{r, f )• In addi- 
tion, the insets demonstrate that Czzij") decays toward 
zero very rapidly, surprisingly independent of the magni- 
tude of spins and angle 9. 

On the other hand, the twisted magnetic fields can di- 
rectly alter spin configurations in the ground state, which 
in principle reflects intrinsic properties in the different 
spin case. We thus analyze how the direction of spin 
polarization changes with site-i by introducing a classic 
polarized angle defined as 



a, = tan 



sn/{s.. 



(4) 



which essentially measures the deviation of the spin po- 
larization at Figure 121 shows with three different values 



a,: 
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-{Tr-9)+Tr mod (i- 1,2) (5) 



where ares is a residual angle which is two orders smaller 
than ai but depends on i, hi and h^r. One can see that ai 
depends almost linearly on i, implying that ground state 
configuration displays a perfect classical spiral structure 
in spite of strong quantum fluctuations. The present 
ground state is just the superposition of some originally 
lowest-lying states whose energy is the order of the edge 
excitation energy and depends also weakly on the mag- 
netic flcld. Since arbitrary boundary fields do not change 
the integrability of the Hamiltonian (1) for 5* = 1/2, there 
may still be spinon-like topological cxcitations|20j. How- 
ever, for 7^ or TT, the total spin and its z— component 
are no longer good quantum numbers, and the spinon 
essentially does not carry definite spin in contrast to the 
9 = case where the spinon possesses a spin of one half. 
This implies that the spinon propagates through some 
kind of spiral path in response to twisted magnetic fields. 
Nevertheless, since the integrability of the system ensures 
the spiral spinon is dissipationless, it is unclear what the 
topological charge of the spiral spinon is and what kind 
of conservation law accounts for the topological charge of 
the spiral spinon, which is beyond the current studies. 
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FIG. 2: The spin polarization angle ai as a function i with 
various twisted angle 6 &t h = 1. In the figure, at are explicitly 
displayed with odd i only, while the results for even i are 
obtained from shifting those for odd i upwards by tt. The 
solid curves on the right panel fit Eq. (6). 

For the S = I case, the polarization angle can be ex- 
pressed with a significantly different form 



ai ~ TT mod {i — 1,2) 
1 



(6) 



-(tt -9)il + sign(^ - -)(1 - exp(-^^))) 



where sign(a;) is the sign function. It is interesting to 
notice that for either odd or even i, ai exhibits a kink at 
i = N/2 with a transient width ^e,iv which in principle 
is proportional to the correlation length for sufficiently 
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large N. However, it depends on the value of and in- 
dependent of the magnitude of external fields as we have 
checked for h = 0.1 and 10 cases. ^g^N = 2.9, 3.6 and 4.0 
for — 7r/3, it/2 and 2tt/3. The kink can be interpreted 
as a consequence of the soliton-like nature of the ma gnon 
excitation owing to the presence of the gap for S' = 1 [lj| . 

For the quantum spin-1 HAF chain, there exists the 
Z2 X Z2 hidden symmetry, which makes its excitations 
unstable against disturbances [23j| . It has been found 
that an open boundary can result in four-fold degeneracy 
with S* = 0, 1 in the thermodynamic limit |22||. while the 
magnetic doping causes a midgap state |2J, |25j. In the 
present case, the boundary field simply lifts the degener- 
acy and induces subsequently a midgap. Figure O shows 
that a midgap energy emerges linearly in h for small h 
and enters the continuum when h > 0.5. Therefore, the 
boundary field can be devised as a tuning parameter to 
manupulate such a macroscopic quantum state , w hich is 
detectable in neutron scattering measurements[2^. 
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FIG. 3: The midgap energy versus h with various 9 for the 
S = 1 case. The horizontal line indicates the Haldane gap. 

Now we turn to the intrinsic transport properties as- 
sociated with above findings. We notice that the twisted 
magnetic fields essentially impose a "spin voltage" be- 
tween two ends of spin chains[ll|. For 5* = 1/2, the spin 
spiral polarization homogenously extends over the whole 
chain, resulting in a spin current which does not depend 
on the chain length (see below). This supports the con- 
juncture of Zotos and coworkers on an ideal metallic spin 
chain with S = 1/2 |^, resulting from the integrabil- 
ity. In this sense, the correlation function Czz{r) in the 
bulk shown in Fig^properly refiects the quasi-long range 
character and the edge effects near the end of the chain. 
However, for the 5* = 1 case, the kink prevents the prop- 
agation of the magnon over characteristic distance 2^e jv- 
One can observe the spin current only when the size is the 
order of 2^6i,Ar or smaller. This is consistent with the con- 
juncture on the persist spin current on a ring proposed 
in Ref . [l^ . Czzir) in Fig. ^indeed deplays a short-range 
correlation for the spin-1 chain, while the tails confined 
nearby r ^ N unveils the fact that edge spins are indeed 
asymptotically free due to the breaking of the hidden 
Z2 X Z2 symmetry. 
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FIG. 4: The spin current multiplied by A'' shows a significantly 
difi'erent size-dependence for both 5 = 1/2 and 1. h equals to 
one for various cases and is explicitly indicated otherwise. 



Following Shen[2] 



we introduce a spin current: 
- (S, X S,+i). 



(7) 



which is site-independent. As the twisted magnetic field 
is applied on the xz plane for both i and N, only y- 
component is non-zero. Figurc^shows NI^ as a function 
of the chain length N for both 5 = 1/2 and 5 = 1. 
When N increases, it increases for small h or unchange 
for large h for S = 1/2, but it always decreases for 5 = 1. 
The results for S = 1/2 indicate again an ideal magnetic 
metal for the spin-half HAF chain, while generally an 
insulator for S = I except for small size systems where 
the spin current is observable |l2j. As far as the spin 
dissipation is concerned in the S = 1 case, we found that 
Is = Ae-^/« with A = 0.39 and ^ = 6.04 apart from a 
weak ^-dependence. ^ is in very good agreement with the 
correlation length of the ordinary S = 1 HAF chain. |2^ 
Since the S = 1 chain is an insulator, we focus on the 
5=1/2 case below. We notice that the spin current is 
very sensitive to the twisted angle 6 as shown in Fig[Sfa) 
and (b). On the other hand, as seen from Fig. EJc), 
the spin current rapidly increases with increasing h and 
reaches a maximum at about h = 1 which is the critical 
value of the uniform field to saturate the uniform mag- 
netization. When h is further increased. Is saturates. 
For very small hi and h^, the system lies in the linear 
response regime. In this case, we found that the spin cur- 
rent is proportional to the vector hi x h2 which precisely 
follows sin 6* as demonstrated in Fig. [Sfb). Therefore, we 
can naturally devise a "spin voltage" in a form of 



Vs = hi X h2 



(8) 



a "spin- vector potential" for Eq. (1). Then we obtain 
the spin current as follows: 



ls=D- V, 



(9) 
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FIG. 5: The spin current and the spin conductance for 
S = 1/2. (a) For various h and (b) For the linear re- 
sponse regime with different /ii and /iat. Solid curves draw 
sin0,which indicate that Is is proportional to hixh2; (c) Is as 
a function of h with 6 — n/2. {d)The same size-dependence 
of the spin conductance with all h as the same as in (b) . 



which is one kind of the London equation with D being 
the spin conductanee. Sinee both sides of the above equa- 
tion are time-reversal invariant so that the spin eurrent 
is dissipationless0. The spin eonductivity characteriz- 
ing the bulk properties described by Eq. (1) is given by 
a = DN. We calculated the spin conductance D for 
various hi and h^r with different chain length N as ex- 
plicitly demonstrated in Fig. I3d) and found that it is 
linearly scaled to zero 1/N. With making extrapolation 
for the thermodynamic limit, we obtain a = 10.0 for 
the Heisenberg antiferromagnetic spin-half chain, which 
is expected in connection with the Drude weight studied 
by Zotos with the finite temperature 0- 

Is the spin current discussed here observable? The sys- 
tem discussed here is implementable experimentally as a 
two leads spin systems introduced recently in Ref . [ll| . A 
spin current generates an electric field. By measuring the 
electric voltage difference between two points in the vicin- 
ity of the spin chain, one can detect the spin current. 
We propose another way of detecting the spin current in 
a spin half chain. With a fixed hi and a rotating h.2, an 
alternating spin current is generated and an alternating 
electric field is subsequently observable nearby the spin 
chain. Measurement of the ac voltage of a given point 
close to the chain and the reference point (ground) re- 
veals the spin current. 

In conclusion, we study the spin transport properties 
of the Heisenberg spin chains for S = 1/2 as well as 
S = I via applying twisted boundary magnetic fields. Al- 
though the boundary conditions generally do not affect 
bulk of a sufficiently large system, the twisted bound- 



ary fields indeed change entirely spin orientation for the 
chains in the thermodynamic limit, allowing to detect 
the spin transport properties. The significantly different 
transport properties are found for S = 1/2 and S = 1 
chains. The former is spin-metallic and has spiral spinon 
excitation, while the later is spin-insulator which involves 
a static kink and unveils a midgap state lift from the de- 
generate ground state by the external field. A London- 
type equation for spin current with the spin voltage and 
the detection of an alternating electric field are proposed 
for the spin-metallic case in the linear response regime. 
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